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Abstract 

The conformational properties of a semi-flexible polymer chain, anchored at one end in a uniform 
force field, are studied in a simple two-dimensional model. Recursion relations are derived for the 
partition function and then iterated numerically. We calculate the angular fluctuations of the 
polymer about the direction of the force field and the average polymer configuration as functions 
of the bending rigidity, chain length, chain orientation at the anchoring point, and field strength. 

PACS numbers: 36.20.Ey, 83.50.-v, 87.15.Aa 
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I. INTRODUCTION 



The influence of external forces on the conformational properties of polymers has been 
studied extensively in recent years. Polymers stretched by attached magnetic beads by 
laser tweezers 0, 0], and by optical fibers ^ and polymers in flow fields |^ have 

received much of the attention |Q. The study of polymer deformation in elongational flow 
goes back to the prediction of a coil-stretch transition |rT| and early birefringence and 
light scattering experiments |T^, r^|. Experimental techniques ||^, T^, W\ which allow 



direct visualization of polymer conformations in simple flows have given this field a new 
perspective. Here the main idea is to use fiuorescently labelled DNA molecules, which are 
long enough so that their conformations can be resolved in an optical microscope. 

In contrast to typical synthetic polymers, DNA chains are semi-flexible, with a persistence 
length of about 80 nm. For contour lengths of a few microns or more these chains behave 
as flexible polymers, in the absence of external forces. In the case of highly stretched DNA 
chains the bending rigidity has been shown to play an important role [l^, |l^ . 

The force-extension curve of a semi-flexible polymer pulled at both ends is derived in 



Ref. [|T^. Predicting the deformation of a polymer in a flow field is considerably more com- 
plicated for two reasons. First, there is a direct hydrodynamic interaction between different 
polymer segments [2^, Second, even if the conformation-dependent, fluctuating 

drag on each bead is approximated by a friction term proportional to the local flow velocity 
("free-draining" approximation), the force on each bead depends on the positions of all other 
beads. Thus, most theoretical studies have rehed on computer simulations R 11^1, 120, pTl 122 



and / or consider flexible chains [T^, |2^, ^ ^ ^ 



In this paper we study the conformational properties of a semi-flexible chain, anchored 
at one end, in two dimensions in a constant force field. In our model the polymer partition 
function is determined by simple recursion relations, which are easily iterated numerically. 
Very little computing time is required, and there is no statistical error in the results, but 



some other approximations, such as the Villain approximation are involved, as will 
become clear. 

The two-dimensional model of a semi-flexible polymer is described in Section II. Recur- 
sion relations for the partition function are derived in Section |T|. In Section |V| we calculate 
the angular fluctuations of the polymer segments about the direction of the applied force, 
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and in Section the longitudinal extension due to the force. In Section |V| we vary the angle 



between the polymer and the force field at the anchoring point and see how this affects the 
mean polymer configuration. Finally, in Section |V11| the case of a polymer pulled at its ends 
is briefly considered. 

II. THE MODEL 

In the wormlike chain model of a semi-fiexible polymer, the Hamiltonian is given by ||2^ 

Ho = ^ £ ds {d,t)\ (1) 

where t is a unit tangent vector and s is the arc length. We consider a discrete version of 
this model in two spatial dimensions, with Hamiltonian 

N 

no = jY,i^i-9i.if. (2) 

i=l 

The polymer chain consists of + 1 line segments of fixed unit length. The ith segment 
forms an angle 6i with the x axis. One end of the polymer is anchored at the origin, and 
the orientation angle 6o of the first segment is also assumed to be fixed. 
To include a uniform force field Fq in the x direction, we add the terms 

N N i N 

Hi = -Fq = ~^<^J2Y1 = -FoJ2i^ + ^-'^) COS ^^ (3) 

i=l i=l j=l i=l 

to the Hamiltonian. The external field could be an electric, gravitational, or uniform flow 
field. The partition function corresponding to Eqs. (0) and (^) is given by 

/CXD poo I ^ I 

dOi... dON exp I- J2\KiO.,-9,_i)^ + L{N + COS Oi \, (4) 

oo J-oo [ i=i -'J 

where K = J/kT and L = F^/kT. 

A nice feature of this model is that it can be solved exactly in the absence of an external 
field, i.e. for L = 0. The mean square end-to-end distance is given by 



1 - e-^/« 

N- 



(5) 



eV? - 1 

with the persistence length ^ = AK. For ^ ^ ^ 1, {R\i) ~ (1 + 2^)A^, corresponding 
to an ideal flexible chain with Kuhn length 1 + 2^. In the limit N <^ ^ with ^ ^ 1, 
{R%) ~ (A^ + 1)^, corresponding to a rigid rod. 



To obtain a more tractable model, we make the Villain approximation p7 



exp 



L cos 9 



^ exp [-L{e - 2nmf] 



(6) 



m=—oo 



in Eq. (^. It was originally introduced in studies of the two-dimensional x — y model and 
the roughening transition, where the Hamiltonians have a similar form. The approximation 
preserves the periodicity of the cosine function but leads to more manageable Gaussian 
integrals. An irrelevant normalization factor on the right-hand-side has been omitted. The 
constant L may be determined by expanding both sides of Eq. (^) in Fourier series and 
equating the lowest two Fourier coefficients. This yields ||2^ 

-1 



4 In 



hCL) 



(7) 



where Jq, /i are modified Bessel functions. 

Replacing X]m=-oo -^1- © Sm=-mniax defines a further approximation, which may 
be systematically improved by increasing mmax- In the finite mmax approximation, con- 
figurations with up to mmax loops about the origin receive the same statistical weight as 
for mjnax = oo, but the statistical weight of configurations with more than mmax loops is 
underestimated . 



III. RECURSION RELATIONS 

As a first approximation we neglect all but the m = term in Eq. (P), replacing e^'^"'^^ by 
g-L6i2 fjj^jg jg good approximation for sufficiently large K and/or L. The corresponding 
partition function is 

poo poo 

Z%ih, ... ,/i^)= / d9^... / d9 



N 

oo J —oo 



N 



X 



exp <^ [-^(^* - ^-i)' - L^^' + \ , (8) 



i=l 



where 



Li = {N + l-t)L, (9) 



as follows from Eq. . The hi are auxiliary variables that will be used in calculating thermal 
averages. 
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The partition function (||) may be evaluated by straightforward integration over 9i,92, ■ ■ ■ ■ 
The first k — 1 integrations contribute 



poo 




/ d9i... 


/ d6k-i exp <^ 


J — oo 


J —oo \ 



where 



lk-1 + K 



exp 



7fc-i + K 



Qk-i 



71 = K + Li, 



Ik 



Ik-i + K 



k = 2,...,N , 



(10) 



= Qk exp [-'jkdl + 2/?fc6'fc] , 
where is a constant, independent of 9k. Equation ( p!0D and the recursive property 

/oo 
d9k-i exp [-K{9k - 9k-iy - 1^91 + hk9k] Ql.M-i) (11) 
-oo 

imply 



(12) 



(13) 
(14) 



/5i = K9o + -/ii 



A; = 2, 



AT 



7/c-i + 2 

Iterating Eq. ([l^) to obtain and using = / d9NQ%{9]^) = 

= (7i"/7Af)^^^ exp [/3^/7Ar] Qn, as follows from Eqs. (||) and (p!0D, we obtain 

1/2 



7° 



TT 



TT 



TC 



X exp 



+ 



+ ■■■ + 



PI 



(15) 
(16) 



(17) 

_77v 77V-1 + ^ 7i + 

The partition function ([T7|) is completely determined by the recursion relations ([T3|)-(p!6D 
and the Li defined in Eq. (pi). 



IV. ANGULAR FLUCTUATIONS 



We now derive the angular fluctuations {9f) of the polymer chain in a constant force field 
from the partition function of Eqs. (^ and (^. In this section we set the initial angle 



^0 and all of the auxiliary variables hi equal to zero. In this case 

K{e^) = f,{L/K,N) (18) 

only depends on K and L in the combination L/K. This follows from rescaling the angles 
9i in the partition function (H). A second consequence is that the Pk in Eqs. (p!5|)-(p!7|) all 
vanish. 

Using the definition (y) of Lj, we write the recursion relation (0) for 7^ in the form 

^„ = K + L{N + l-k)~ . (19) 

According to Eqs. (|^), (0), and (|19D the angular fluctuations {6f) = — din Z%/dLi satisfy 

'«> = 77^ ■ (21) 



32 ^ 
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We have calculated {9^) by numerical iteration of Eqs. (P^, (|T^, (0), and (PT]). The 
results for L/K = 0.1 and = 100, 500, 1000 are shown in Fig. 0. For these three values of 
A^ the results for {6]^) in Fig. |I]are practically indistinguishable. The same is true of {6'jf_-^), 
(^Af-2)) ^^c- The values A^ = 100, 500, and 1000 are large enough so that the angular 
fluctuations at the free end of the chain are independent of the chain length. 

We now examine the A^-dependence of the fluctuations. In Fig. ^ the quantity K{6%) is 
plotted as a function of A^ for five different values of L/K. There is an obvious crossover 
from A^-dependent to A^-independent behavior as A^ increases. We denote the approximate 
value of A^ at the crossover by A'min- According to Eq. (p!8|), A'min only depends on K and L 
in the combination L/K. Fig. ^ shows Amin as a function of L/K. The data are in excellent 
agreement with 

AT^in ~ {L/K)-'/\ L/K « 1 . (22) 

The power law (^) follows from the following argument: For L/K <^ 1 the recursion 
relation (^) implies 

K 

K 

where satisfies the recursion relation 



l + kA,^ + 0{{L/Kf) 



(23) 



Ak={'-^^ Ak^i + {N + l-k) (24) 
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with Ai = N. Since the factor multiplying Ak-i in Eq. (p^ approaches unity as k increases, 
An ^ N"^ for ^ 1. For fixed L/K <^ 1 and small, the first term in 77V, as given 
by Eq. (p3D, is clearly the dominant term, and 7Ar — K/N. However, the second term 
becomes increasingly important as increases. It is reasonable to assume that a crossover 
to a different L/K dependence occurs when the second term in 7Ar, as given by Eq. (P^D, 
becomes comparable with the first, i.e. NA^ ~ A^^ ~ K/L for A^ ~ A^mm- This leads to 
Eq. (H) and the prediction 7^ ~ 77v^.^ ~ K{L/Kf'^ for L/K < 1, A^ > A/'^in- 
In the complementary regime L/K ^ 1 the recursion relation (|19D implies 

L K 

K^'-2L 

for arbitrary A^. Since this result is entirely independent of A^, the large A^ behavior has its 
onset at 



In 



K 



+ O {{K/Lf) 



(25) 



AT^in ~ 1, L/K > 1 . 



(26) 



We now derive exact analytic expressions for the fiuctuations (6*^), ... at the end 

of an infinitely long chain. Reverse iteration of Eq. (|T^, beginning with 7^?, leads to the 
continued fraction 



jN = K + L 



K' 



K' 



K' 



K^ 



(27) 



2K + 2L- 2K + 3L- 2K + AL - 2K + NL 

In the limit N ^ 00 the continued fraction may be evaluated with the help of Eq. (9.1.73) 
in Ref. ||2^, yielding 

-1 



7o 



K 



- 1 



2K 



V 



(28) 



Here J^l,{^z) is the standard Bessel function, and J'J^z) is its derivative with respect to z. 



From Eq. (9.3.23) in Ref. and Eq. (|3), we obtain 



7c 



i/3r(2/3) (L 



3I/3 



Kl 



r(i/3) 



1/3 



L 
K 



< 1 



(29) 



— 1- — 
K K^ 2L' 

These limiting forms are consistent with the expressions for 700 ~ 7Afmin small and large 
L/K given in Eq. (^) and the paragraph which preceeds it. 



/^From the result ( pHD for (6'^) = (27Ar) ^ in the large- limit, it is straightforward to 
calculate (^ir^i), (^^_2), • • • using Eqs. (|19|)-(|1]). 

In Fig. |, K{e%) = K/{2^n) is plotted as a function of L/K for N = 10, 100, 1000, and 
compared with the analytic prediction ( [28| ) for N ^ oo. 

As stated in Section |T|, the mmax = approximation is accurate for sufficiently large K 
and/or L. One can use the results of this section to determine the domain of validity more 
precisely. The approximation, i.e. replacing cos6' by 1 — ^6'^, should be quite reliable if, say, 
{^n) = (27Ar) < (7r/4)^ or > S/vr^ ~ 1. Computing •jn by numerical iteration, one 
can readily check whether this inequality is satisfied for particular values of K, L, and A^. 
According to Eqs. (^), (^1), and (|25D , the inequality 7Ar > 1 corresponds to K > N for 
L/K < 1 and A^ < N^i^, to K{L/Kf/^ > 1.4 for L/K < 1 and A^ > A^min, and to L > 1 
for L/K ^ 1 and arbitrary A^. 



V. LONGITUDINAL STRETCHING 



For a polymer in a constant force field, one of the main quantities of interest is the average 



extension in the fiow direction. According to a prediction of Marko and Siggia [|1S 



^ 1 - CiLKNy^' , (30) 

with Co = 1. This result was derived by approximating the restoring force at position s 
along the chain with the thermodynamic result for a polymer of length s pulled at its ends. 
Note, however, that a polymer in a fiow field fiuctuates most strongly at the free end and 
not at all at the anchored end, quite unlike a polymer pulled at its ends. The derivation also 
assumes a boundary condition (cos On) ^ LK"^ ^ 1 at the end of the chain, at odds with the 
exact result {6%) ~ [LK'^)^^^^ in Eqs. ( ^I]) and (|29|) for large K and A^, where our discrete 



model is equivalent to the continuum model of Ref. |T8[. One advantage of our approach 
is that it avoids these assumptions and yields numerically exact results for the model with 
partition function (^). 

We have checked Eq. (^OD for our model, using the relations 

N N 

(xn) = $^(cos^^,) = - ^(e^^^ + e-^^) , (31) 
i=o j=o 

— -^Af(0; ■ ■ ■ ) hj = ±i, 0, . . . ) . , 

^ ~ ^^(o,...,o) • ^^'^ 



Here the numerator on the right side of Eq. (^) is the partition function of Eqs. (H) and 
(p!7|) with all of the auxiliary fields set equal to zero except that hj = i. In the denominator 
all of the auxiliary fields, including hj, vanish. 

Calculating these partition functions recursively using Eqs. (|T3|), (0), and (|r^), with the 
Li defined by Eq. {^), we obtained the results for the extension shown in Fig. |^. The data 
for sufficiently large K do indeed confirm Eq. (jS^). For our model Co — 0.23. 

It is instructive to compare the chain length with Nmin in the regime where Eq. (|30[) 
applies. For K = 1000 and = 10 000 (filled circular points in Fig. |]) the six-decade 
interval 10° < LKN < 10^ corresponds to 10"^° < L/K < 10'^, or (see Eq. (H) and Fig. D 
to 20 000 > Nmin > 200. Thus, the inequality ^ A^mm only holds for the last three 
decades of LKN. 

There are deviations from the straight line in Fig. |^ for smaller K and larger L. For 
L^ K, the recursion relation ([T9|) implies 



7fc^(Ar + l-fc)L, (33) 

and 



AT ^ N , TV 



{x^) = J2i^ose,) ^ AT - 1 ^(^5 ^ ^ - ^ E - • (34) 

i=0 i=0 1=0 

The sum over the chain segments is easily carried out and yields 

InA^ 

M ^N-^. (35) 

In Fig. ^ numerical data for several bending rigidities, chain lengths and force fields are 
compared with this result. The agreement is excellent. Thus, for a sufficiently strong force 
field or a sufficiently small bending rigidity, the chain extension varies as L~^, just as for 
flexible chains with arbitrary L. 



VI. MEAN POLYMER CONFIGURATION FOR 6*0 7^ 

In this section we consider polymer conformations with the first segment fixed at a non- 
zero angle 6q with respect to the direction of the force field. We make the Villain approxima- 
tion and restrict the sum in Eq. (|^) to m = —1, 0, 1. In the corresponding polymer partition 
function there is a separate sum over m for each polymer segment. For a sufficiently stiff 
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polymer and/or a sufficiently strong force field, the conformations of the chain are domi- 
nated by the same potential minimum. In this case the separate sums may be replaced by 
a single sum, leading to the simpler partition function 



/oo roo 
dOi... I dOt 

X exp jj^ [-K{e, - 9,_,f - 1,(9, - 2nmf + JiA] | • (36) 

Here we have again introduced a set of auxiliary variables hi, to be used in constructing 
thermal averages. 

Expanding {6i — 27rm)^ in powers of 6i, one sees that the partition function ( |36D can be 
expressed as 

1 

ZN{hi, .■.,~hN)= Yl e-^"''^'(^i+-+^'^)Z^(/ii + 47rmLi, . . . ,hN + ^nmLN) (37) 



m=—l 



in terms of the m = partition function Z'^lhi, . . . , Hn) defined in Eq. (|T^. Since we know 



how to calculate numerically with the recursion relations of Section |nT|, we can also 



calculate via Eq. (pTf). Using Eqs. (H), ([T3D -(|17|), (^), and the relations 

k k 



{x,) = J2(^ose,) = i X^(e^^^ + e-^'^) , (38) 

3=0 j=0 

(y,) = X^(sin0,) = i J2^e^'' ' ^"'') ' (^9) 
i=o i=o 

/g±ie,\ ^ ZN{0,...,hj = ±i,0,...) 

^ ^iv(0,...,0) ' ^ ^ 

we have evaluated the average position of the jth segment of the polymer chain for fixed ^o- 

Figs. show how the average position depends on the tilt angle 6^0, field strength L, 
and bending rigidity K. In all of these figures = 100. 

Figs. and ^ show the average polymer configuration in the x, y plane. As the field 
strength increases, the polymer is bent towards the field direction and is stretched longitu- 
dinally. The elongation is more pronounced for smaller bending rigidities. 

In Fig. P the transverse extension of the polymer as a function of the tilt angle is shown 
in more detail. The curves are sinusoidal for small field strengths L but for larger L bend 
abruptly near 6*0 = ivr, due to the instability of a polymer directed against the force field. 
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Since our model includes fluctuations, the polymer "tunnels" between the two equivalent 
free-energy minima, and there is no spontaneous symmetry breaking ai 6q = n. 



Fig. |T0| shows the tranverse extension as a function of the field strength for three different 
Oq. For = 10, = 100, and L < 10^^, the effect of the force field is negligible. For 
stronger force fields, there seems to be a regime where {yj\i) ~ L^^/^. 

The contour length Ci of the average configuration (see Figs. |^ and is shown in Fig. |TI|. 
Again, the effect of the force field is negligible for = 10, = 100, and L < 10~^ . Varying 
the tilt angle only affects the contour length near the onset of the deformation due to the 
force field. 

Finally, we have considered the angle 6 between a line through the end-points of the chain 
and the direction of the force field. A weak force field deforms the polymer only slightly, 
and 6 varies linearly with the tilt angle Oq. For a strong force field, on the other hand, 6 
varies abruptly as 6q approaches vr, due to the instability mentioned above. The behavior 



of (5 as a function of L in Fig. O is qualitatively similar to that of {un) in Fig. 10. 



VII. POLYMER PULLED AT ITS ENDS 



Thus far we have considered an external force field that acts on each monomer of the 
semi-fiexible polymer. With only minor modifications the case of a constant force applied 
at the ends of the polymer can also be studied. Anchoring one end at the origin at a 
fixed angle 6*0, we apply a constant force at the other end by replacing Eq. (||) with Hi = 
—FqXn = — -fo X^j^i "^os ^j- Villain approximation the partition function is again 

given by Eq. (||), but Eq. @ is replaced by Li = L. With this definition of the Li the 
partition function may be calculated recursively, as in Section |IT1[ For symmetric boundary 
conditions at the ends of the chain, i.e. 6o and 6n both fixed or both free to fluctuate, the 
calculation is also straightforward. 

For fixed = and fluctuating 6]\f, the angular fluctuations of the polymer segments are 
given by Eqs. (|TB|), (^, and (^, with Eq. replaced by 

j, = K + L--— . (41) 

K + Ik-i 
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In the long chain hmit approaches the fixed point 

7oo = §+ (^] +KL (42) 



of Eq. flD. 

For large K and our discrete model is equivalent to the continuum model of Marko 
and Siggia [0. In this limit Eqs. (plf ) and imply the same result {6%) = {AKL)^^/'^ for 



the angular fiuctuations as in Ref. |T8[ . 

In Fig. 0, K{e%) = K/{2-fN) is plotted as a function of L/K for = 10, 100, 1000, and 



10 000 and compared with the analytic prediction (|4^) for ^ oo. 

We have also calculated the mean configuration of a polymer pulled at its ends for fixed 



^0 > and fiuctuating 9^. Fig. |IJ shows the tranverse extension as a function of the force for 
three different Oq. For K = 10, N = 100, and L < 10^^, the effect of the force is negligible. 
For stronger forces there seems to be a regime where (yiy) ~ We found quite similar 

behavior for a polymer in a uniform force field, as shown in Fig. 0. 

VIII. CONCLUDING REMARKS 

For calculating the conformational properties of a semi-fiexible chain in a uniform force 
field our recursive approach has several advantages: (i) It requires very little computing 
time, and (ii) it allows one to consider very long chains. For a clearly-defined model exact 
numerical results are obtained. Thus, (iii) there is no statistical error, and (iv) some of the 
approximations in earlier theoretical work are avoided. Finally, (v) the recursion relations 
furnish some analytical insight. We were able to obtain some exact results for the asymptotic 
behavior of long chains. While most previous studies have focused on the force-extension 
relation, we have also analyzed angular fiuctuations. 

A disadvantage of the approach is the limitation to two dimensions. However, many 
of the results probably apply, at least qualitatively, to chains in three spatial dimensions. 
Furthermore, the results are directly applicable to polymers confined to two dimensions, for 



example, DNA electrostatically bound to fiuid lipid membranes |]3^ . 

The Villain approximation was used to obtain a tractable model. It preserves the peri- 
odicity in 6 and is no more unrealistic than using a quadratic bending energy for arbitrary 



12 



angles or ignoring excluded volume. We only presented results for the single-m approxi- 
mation with mmax < 1, which underestimates the statistical weight of configurations with 
segments pointing in widely different directions but is accurate for sufficiently large K and/or 
L. 

The single-m approximation can, of course, be improved at the cost of greater computing 
time. Retaining all the Villain sums leads to the partition function 

oo oo 
mi=— oo mjv=— oo 

X Zli{hi + A-KTYiiLi, . . . , /iTv + 47rm7vLAr) (43) 

instead of Eq. (P7|). Here is the m = partition function in Eq. (|T^, which we know how 
to compute recursively. Usually one is interested in K and L for which the angular differences 
between adjacent segments are small, certainly less than 27r. Then the sums on the right 
side of Eq. ( ^3|) may be restricted to the terms with — 1 < mi < 1, nii — 1 < m2 < mi + 1, 
etc. Computing with no further approximations requires 3^ evaluations of Z^. 
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Figure Captions 



Figure 1: Average angular fluctuations K{9f) for L/K = 0.1 and N = 100 (•), = 500 
W, = 1000 (*). 

Figure 2: Average angular fluctuations K{9%) of the final chain segment as a function of A^ 
for L/K = 10-5 (•), L/K = 10"^ (o), L/K = IQ-'^ (A), L/K = IQ-^ L/K = 10"^ 
(*). The straight line has slope 1. 

Figure 3: A^^in as a function of L/K for AT = 1 (•), AT = 10 (o), AT = 100 (A), AT = 1000 
(^). The straight line has slope -1/3. 

Figure 4: K{9'j^) = K/{2'yj\f) for a polymer in a constant force field as a function of L/K 
for A^ = 10 (•), A^ = 100 (o), and A^ = 1000 (A), together with the exact result (H) 
for N oo {-k). The straight lines have slopes —1/3 and —1, respectively. 

Figure 5: Average extension, (xjv) of the chain in the direction of the force field direction 
as a function of the field strength A, for AT = 10, A^ = 100 (*); K = 10, A^ = 1000 
(o); K = 10, N = 10000 (o); K = 100, A^ = 1000 (A); K = 100, A^ = 10000 
K = 1000, A^ = 10 000 (•). The straight line has slope -1/2. 

Figure 6: Average extension (x^) of the chain in the direction of the force field as a function 
of the field strength L. The symbols correspond to the same parameters as in Fig. |^. 
The straight line has slope —1. 

Figure 7: Average positions (xj), of a chain of length A^ = 100 with A' = 1, = ^/3 
(lower panel), and 6o = 2tt/3 (upper panel), with L = 10^^ (o), L = 10~^ (A), 
L = 10~^ {-k), L = 10"^ (•). The straight line indicates the tilt angle Oq. 

Figure 8: Average positions (xi), {yi) of a chain of length A^ = 100 with A' = 10, = 7r/3 
(lower panel), and 9q = 2n/3 (upper panel), with L = 10~^ (o), L = 10~^ (A), 
L = 10~^ {-k), L = 10^^ (•). The straight line indicates the tilt angle Oq. 

Figure 9: Average transverse extension (yjv) as a function of Oq for A^ = 100, A' = 10 and 
L = 10-9 (*), L = 10-5 (*) , and A = 10-^ (A) , A = 10-^ (o) , A = 10-^ (•). 
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Figure 10: Average transverse extension (t/at) as a function of L for = 100, K = 10, and 
9o = n/3 (•), 00 = 2n/3 (o) , 9^ = ITvr/lS (A). 

Figure 11: Contour length Ci of the average configuration as a function of L for = 100, 
K = 10, and Oq = n/lSO (•), Oq = n/3 (o), Oq = 2n/3 (A). 

Figure 12: Angle 6 between a hne through the end points of the chain and the direction of 
the force field as a function of L for A^ = 100, K = 10, and Oq = vr/lSO {•), 6o = tt/3 
(o), ^0 = 271/3 (A). 

Figure 13: K{9'^) = K/{2'~^n) for a polymer pulled at the ends as a function of L/K for 
A^ = 10 (*), A^ = 100 (A), A^ = 1000 (o), and A^ = 10 000 (•), together with the exact 
result (^2]) for A^ ^ oo (^). The straight lines have slopes —1/2 and —1, respectively. 

Figure 14: Average transverse extension (y^v) of a polymer pulled at its ends as a function 
of L for A^ = 100, K = 10, and Oq = tx/3 (•), = 2tx/3 (o), Oq = ITvr/lS (A). 
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FIG. 1: Average angular fluctuations K{ef) for L/K = 0.1 and N = 100 (•), N = 500 (*), 
N = 1000 (*). 
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FIG. 2: Average angular fluctuations K{9ff) of the final chain segment as a function of N for 
L/K = 10-5 (•), L/K = 10-^ (o), L/K = 10"^ (A), L/K = 10"^ L/K = 10"^ (*). The 
straight line has slope 1. 
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FIG. 3: N^in as a function of L/K for K = 1 {•), K = 10 (o), K = 100 {A), K = 1000 (★). The 
straight hne has slope -1/3. 
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FIG. 4: K{9'j^) = K/{2^]\f) for a polymer in a constant force field as a function of L/K for = 10 
(•), N = 100 (o), and N = 1000 (A), together with the exact result (|8|) for TV ^ oo (*). The 
straight lines have slopes —1/3 and —1, respectively. 
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FIG. 5: Average extension, [xm) of the chain in the direction of the force field direction as a 
function of the field strength L, for K = 10, N = 100 (*); K = 10, N = 1000 {o); K = 10, 
AT = 10000 (o); K = 100, N = 1000 (A); K = 100, N = 10000 (*); K = 1000, N = 10 000 (•). 
The straight line has slope —1/2. 
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FIG. 6: Average extension {xn) of the chain in the direction of the force field as a function of the 
field strength L. The symbols correspond to the same parameters as in Fig. |5|. The straight line 
has slope —1. 
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FIG. 7: Average positions (xj), (j/j) of a chain of length TV = 100 with K = 1, Oq = tt/3 (lower 
panel), and 6*0 = 2tt/3 (upper panel), with L = 10"^ {°), L = 10"^ (A), L = 10"^ L = 10"^ 
(•). The straight line indicates the tilt angle ^o- 
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FIG. 8: Average positions (xj), (yj) of a chain of length N = 100 with K = 10, = 7i"/3 (lower 
panel), and 9o = 27r/3 (upper panel), with L = 10'^ (o), L = 10"^ (A), L = 10"^ (★), L = 10"^ 
(•). The straight line indicates the tilt angle ^o- 
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FIG. 9: Average transverse extension (un) as a function of for N = 100, K = 10 and L = 10 ^ 
(*), L = 10-5 (★) , and L = 10"^ (A) , L = 10"^ (o) , L = 10-^ (•). 
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FIG. 10: Average transverse extension (yjv) as a function of L for TV = 100, K = 10, and Oq = ir/S 
(•), 00 = 27r/3 (o) , eo = 177r/18 (A). 
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FIG. 11: Contour length C/ of the average configuration as a function of L for AT = 100, K = 10, 
and ^0 = 7r/180 (•), = 7r/3 (o), 9o = 2tt/3 (A). 
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FIG. 12: Angle S between a line through the end points of the chain and the direction of the force 
field as a function of L for TV = 100, K = 10, and = 7r/180 (•), = tt/3 (o), = 2tt/3 (A). 
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FIG. 13: K{9'j^) = K/{2^j\f) for a polymer pulled at the ends as a function of L/K for = 10 ( 
N = 100 (A), N = 1000 (o), and N = 10 000 (•), together with the exact result (||) for N 
(*). The straight hues have slopes —1/2 and —1, respectively. 
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FIG. 14: Average transverse extension (yj\r) of a polymer pulled at its ends as a function of L for 
= 100, K = 10, and = tt/S (•), = 27r/3 (o), Oq = 177r/18 (A). 
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